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LONG  TERM  GOALS 

The  long  term  goal  of  this  project  is  to  develop  efficient  inversion  algorithms  for  successful 
geoacoustic  parameter  estimation,  inversion  for  sound-speed  in  the  water-column,  and  source 
localization,  exploiting  (fully  or  partially)  the  physics  of  the  propagation  medium.  Algorithms  are 
designed  for  inversion  via  the  extraction  of  features  of  the  acoustic  field  and  optimization.  The 
potential  of  analytic  approaches  is  also  investigated. 

OBJECTIVES 

•  Achieve  accurate  and  computationally  efficient  inversion  for  propagation  medium  parameters  and 
source  localization  by  designing  estimation  schemes  that  combine  acoustic  field  and  statistical 
modeling. 

•  Develop  methods  for  passive  localization  and  inversion  of  environmental  parameters  that  select 
features  of  propagation  that  are  essential  to  model  for  accurate  inversion. 

•  Implement  Bayesian  filtering  methods  that  provide  dynamic  and  efficient  solutions  for  the  first  two 
objectives. 

•  Develop  analytic  techniques  for  sediment  sound  speed  estimation. 

APPROACH 

We  first  worked  on  a  new  sediment  sound  speed  estimation  scheme  based  on  Stickler’s  inverse 
problem  approach  [1,2].  As  we  will  mention  in  the  next  section,  we  managed  to  estimate  sediment 
sound  speed  and  thickness  in  a  synthetic  waveguide  using  an  analytic  method,  improving  on 
previously  obtained  results  with  a  similar  but  simpler  approach. 


1 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

30  SEP  2014  2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2014  to  00-00-2014 

4.  TITLE  AND  SUBTITLE 

Efficient  Inversion  in  Underwater  Acoustics  with  Analytic,  Iterative,  and 
Sequential  Bayesian  Methods 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

New  Jersey  Institute  of  Technology, Department  of  Mathematical 

Sciences, Newark, NJ, 07102 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

_ _ _  ABSTRACT 

18.  NUMBER  19a.  NAME  OF 

OF  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Same  OS 

unclassified  unclassified  unclassified  Report  (SAR) 

7 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


Also,  continuing  efforts  from  previous  years,  we  worked  with  Bayesian  approaches  applied  to  sound 
signals  for  the  extraction  of  acoustic  features  using  a  combination  of  physics  and  statistical  signal 
processing.  We  extended  work  on  the  development  of  a  quasilinear  model  for  source  location, 
bathymetry,  and  water  column  sound  speed  profile  estimation.  We  demonstrated  via  a  comparison  to 
simulated  annealing  that  our  approach  is  superior  to  global  optimization:  although  the  method  is  much 
faster,  the  accuracy  in  the  estimation  is  not  compromised. 

More  work  with  a  new  statistical  model  more  accurately  describing  dispersion  in  the  ocean  was 
conducted  for  identifying  arrival  times  and  amplitudes  of  distinct  frequencies  (within  a  single  mode  or 
across  different  modes).  Those  characteristics  provide  significant  information  on  properties  of  the 
propagation  medium  and  source  location.  The  importance  of  modal  arrival  times  and  amplitudes  in 
geoacoustic  inversion  and  source  localization  using  dispersion  curves  has  been  extensively  discussed 
in  [3,4,5]. 

WORK  COMPLETED 

Sediment  sound  speed  was  estimated  using  a  new  efficient  analytic  method.  Despite  low  frequencies 
used  for  the  inversion,  high  resolution  was  attained  with  sediment  sound  speeds  and  thicknesses 
accurately  extracted. 

We  extended  work  on  the  development  of  our  quasilinear  model  for  source  location,  bathymetry,  and 
water  column  sound  speed  profile  estimation  and  demonstrated  that  it  is  superior  in  terms  of  efficiency 
to  global  optimization  in  tenns  of  efficiency  by  comparing  it  to  fast  simulated  annealing. 

More  work  was  conducted  on  extracting  dispersion  curves  for  long-range  propagation  with  a  new 
statistical  model.  Results  show  that  the  new  model  offers  higher  reliability  in  the  estimation  of 
dispersion  than  conventional  spectrogram  methods  and  allows  us  to  perform  an  accurate  inversion  for 
sediment  sound  speed.  A  new  and  more  accurate  statistical  model  than  the  simple  Gaussian  one  that 
had  been  previously  used  was  developed. 

RESULTS 

Figure  1  shows  sound  speed  estimation  in  sediments  and  the  corresponding  sediment  thickness 
employing  the  method  mentioned  in  the  past  section.  The  interpolant  approach  developed  in  our  work 
provides  closer  profiles  to  the  true  one  than  other  approaches,  with  the  improvement  mostly  evident  at 
the  interfaces  between  different  sediments.  This  feature  is  the  one  that  allows  accurate  estimation  of 
the  sediment  thickness. 
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Figure  1:  Sediment  sound  speed  and  thickness  estimation  with  an  analytic  method,  improving  on 
results  previously  obtained  from  similar  techniques  with  two  different  profiles.  Profde 
discontinuities  are  well  tracked,  providing  estimates  both  for  sound  speed  and  sediment  thickness. 

The  same  method  was  applied  to  noisy  data.  Several  numerical  issues  were  addressed  with  the 
development  of  a  suitable  regularization  technique.  The  final  results  show  that,  even  in  noisy 
environments,  the  new  method  can  estimate  sediment  sound  speed.  Sediment  thickness  was  also 
estimated  in  most  cases.  In  other  cases,  smoothing  occurred,  blurring  the  boundaries  between 
sediments.  Results  are  presented  in  Figure  2  for  three  different  Signal-to-Noise  Ratios.  Experiments 
were  performed  with  a  variety  of  environments  with  similar  results  [6]. 
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Figure  2:  Performance  of  the  analytic  inversion  method  at  three  different  noise  levels. 
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Using  Shallow  Water  06  data,  we  compared  our  linearization/particle  filtering  approach  to  inversion 
using  simulated  annealing.  A  comparison  for  two  sets  of  arrival  times  for  three  paths  (direct,  first 
surface  reflection,  and  first  bottom  reflection)  is  presented  in  Table  1.  The  results  are  practically 
identical.  However,  our  method  required  six  iterations  per  case  whereas  the  simulated  annealing 
method  required  from  a  few  hundred  to  a  few  thousand  iterations  [7].  Figure  3  illustrates  the  simulated 
annealing  estimation  process,  showing  the  number  of  iterations  necessary  for  convergence.  The  new 
feature  of  including  EOF  coefficients  in  the  linearization  process  facilitated  the  estimation  of  water 
column  sound  speed  profiles. 

Table  1:  Comparison  of  the  particle  filtering/linearization  method  and  simulated  annealing  (range 
r,  source  depth  zs,  water  column  depth  WD,  EOF  coefficients  pi,  p.2,  and  p$,  and  tilt. 


Parameter 

Lin.  Estimate 

SA  Estimate 

Lin.  Estimate 

SA  Estimate 

r  (m) 

228 

232 

229 

228 

zs  (m) 

26 

26 

27 

27 

WD  (m) 

76 

76 

77 

77 

Ai 

-54 

-56 

-42 

-43 

^2 

-10 

-12 

-7 

-7 

U3 

2 

1 

2 

-2 

tilt  (deg) 

0 

0 

0 

0 
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Figure  3:  Simulated  annealing  results  for  (a)  and  (b)  source  location,  (c)  bathymetry,  and  (d),  (e), 

and  (f)  EOF  coefficients. 


For  the  dispersion  problem,  dispersion  curves  were  extracted  from  broadband  data,  with  dispersion 
tracks  identified  even  for  weak  modes,  as  shown  in  Figure  4. 
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Figure  4:  Spectrogram  of  a  noisy  synthetic  signal  and  the  dispersion  tracks  as  estimated  by  a 

particle  filter  superimposed. 

We  then  propagated  probability  density  functions  (PDFs)  of  arriving  frequencies  at  particular  times  for 
two  modes  through  a  nonnal  modes  model  to  identify  which  sediment  sound  speeds  would  generate 
the  measured  time  differences;  that  is,  we  exploited  group  velocities  for  inversion.  At  the  output  of  our 
method,  we  obtain  a  PDF  of  the  sediment  sound  speed.  Results  are  shown  in  Figure  5.  The  impact  of 
noise  and  mode  misidcntification  on  the  results  was  also  investigated  [8]. 
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Figure  5:  (a)  and  (b):  PDFs  of  two  arriving  modal  frequencies  at  a  particular  time  as  estimated 
by  the  Particle  Filter,  (c)  The  PDF  of  the  sediment  sound  speed  using  the  frequency  PDFs 

of  (a)  and  (b). 


IMPACT 

The  development  of  the  new  inversion  has  a  significant  impact  on  the  problem  of  geoacoustic 
parameter  estimation,  because  of  its  potential  to  estimate  sound  speed  in  sediments  in  an  exact,  fast, 
analytic  fashion.  The  method  relies  on  limited  prior  information  (sensitivity  to  assumptions  is 
currently  being  further  investigated). 
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Also  the  importance  of  dispersion  analysis  in  geoacosutic  inversion  has  been  long  established,  with 
excellent  results  demonstrated  with  real  data.  Our  method  appears  to  produce  superior  results  in 
comparison  to  conventional  methods  and  provides  narrow  PDFs  of  sediment  sound  speed. 

Lastly,  our  inversion  method  with  multipath  arrival  times  and  linearization  is  fast  and  accurate, 
allowing  for  effective  inversion  in  real  time. 

RELATED  PROJECTS 

The  PI  is  collaborating  with  Drs.  Yardim  and  Gerstoft  on  sequential  filtering  in  ocean  acoustics.  The 
PI  is  also  collaborating  with  Dr.  Leon  Cohen  on  comparing  numerical  and  analytical  descriptions  of 
dispersion  in  ocean  environments  and  is  involved  in  discussions  with  Dr.  Ross  Chapman  on  inversion 
for  attenuation. 
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